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O ■ l INTRODUCTION 

Chandra's excellent imaging resolution allows us an un- 
, precedented glimpse of the spatially resolved X-ray emis- 
5^ ■ sion of nearby galaxies. Unlike the simple nuclear point- 
source emission characterizing the X-ray view of luminous 
quasars, nearby starburst galaxies and active galactic nu- 
clei (AGNs) show a plethora of diverse X-ray features: off- 
center point source emission from X-ray binaries, emission 
from supernova remnants, diffuse collisionally ionized emis- 
sion from galactic halos, emission from photoionized narrow 
line regions (NLRs), as well as the X-ray emission of jets 
and outflows, inc l uding starburst - driven superwinds (e.g.. 



ABSTRACT 

We study the extended nuclear emission of the starburst galaxy NGC 1365. A weak 
obscured AGN and a strong starburst both contribute to the observed X-ray, optical, 
infrared, and radio emission in the inner 2 kpc. The X-ray emission is spatially re- 
solved, allowing comparison with multiwavelength data that highlights the structures 
dominating the nuclear region: the AGN, the nuclear spiral, the circumnuclear star- 
burst ring, and nuclear outflow. The ultrasoft X-ray emission < 0.5 keV is spatially 
coincident with the conical outflow traced by higher excitation optical emission lines 
like [O III] and [Ne III]. The strong starburst concentrated in super-star clusters in 
a circumnuclear ring with radius wlkpc dominates the 0.5-1.5 keV emission and is 
visible in radio, molecular CO, and infrared maps of the central kiloparsec. The hard 
(2-10 keV) emission is dominated by the obscured AGN, but also contributes to the 
emission from relatively old (~ 7My r) but still enshroud ed in dust and extremely 
massive (10 7 M Q ) super-star clusters (|Galliano et al.l 2008). hidden from view in the 
optical and soft X-ray bands. In the Appendix we present the X-ray spectroscopy and 
photometry of BL Lac MS 0331.3— 3629, a high-energy peaked BL Lac candidate at 
z = 0.308, serendipitously detected in one Chandra and five XMM-Newton observa- 
tions of NGC 1365. 

Key words: galaxies: individual: NGC 1365 - galaxies: starburst - X-rays: galaxies. 
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there are at least three 
The obscured act iv e nu- 
iKomossa fc Schulzl Il998l : 
the soft, extende d emis- 
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with optical, radio, and infrared emission maps, the spatially 
resolved X-ray emission allows us to identify important pro- 
cesses associated with the birth and deaths of stars, the size, 
place, and composition of regions of recent star formation, 
as well as to gauge the influence of the AGN on the nuclear 
regions of galaxies. 

The barred supergiantQ spiral NGC 1365 has been 
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sec -2 isophotal radius of NGC 1365 is 30 kpc. 



NGC 1365 is special in the sense that several of 
its multi-wavelength properties are, at first glance, un- 
usual. (1) The degree of excitation of narrow emission- 
lines detected in optical spectra is actually lower in the 
core region of NGC 1365 than in the extranuclear regions 
(e.g., lEdmunds et~al]| 19881 : ISchulz et al II 19991 : IVeilleux" et al.l 
2003). (2) While this could be explained by strong core ex- 
tinction, a bro ad component nev e rtheless shows up i n the 
Balmer lines llVeron et all Il98d : lEdmunds fc Pagell Il982l : 
ISchulz et a.1 .1 1 1 9991 ^1 . The width of the broad-line component 
is ~1600-1900 kms -1 and would formally imply a classifica- 
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Figure 1. Left: Optical VLT image of NGC 1365 with the inner spiral region marked. Top right: HST WFPC2 F555W image of the 
inner spiral of NGC 1365 with the 0.3-10 keV band logarithmically-spaced contours from the 2002 Chandra observation overlaid in red. 
Bottom right: The circumnuclear starburs t ring of NGC 1365 is shown here over the same HST i mage, traced by the optical super-star 
clusters (the SSCs are denoted by sq uares ; [Kristen etalj 19971 . KJL97), the radio hotsplots (circles; Stevens ct al. 1999, SFN99), and the 
12 CO molecular hotspots (ellipses; [Sakamoto et al.ll2007l SHMMP07). 



tion as narrow-line Seyfert 1 galaxy, even though the other 
properties of NGC 1365 do not match that classi fication. (3) 
Com pared to other samples of Seyfert galaxies dWard et al.l 
1988), the optical and IR properties of NGC 1365 appear 
to be off the correlations and tren ds in multi-wavelength 
diagrams (jKomossa fc Schulj ll99St ). (4) To this adds the 
recent finding that the X-ray obscuration of NGC 1365, 
which is unusually variable on very short timescales (e.g., 
iRisaliti et al] 2005a), is likely c aused by the broad- line re- 
gion clouds (|Risaliti et~ai]|2008l) . 

With the key goals in mind of spatially resolving the 
soft X-ray emission seen with ROSAT and disentangling 
AGN and starburst components, getting new clues on the 
multi- wavelength properties of NGC 1365, and studying its 
population of X-ray sources including the luminous source 
NGC 1365-X1, we obtained a Chandra ACIS-S observation 
of NGC 1365 on MPE guaranteed time. 

In this paper we report on the extended X-ray emis- 
sion of NGC 1365. We show that the 0.5-10 keV ex- 
tended nuclear emission originates in the star-forming 



ring and nuclear spiral in the inner 2kpc of NGC 1365, 
as traced by radio hotspots, molecular CO emission, 
mid-infrared (MIR) emission, and optically detected su- 
per star clusters, while u ltrasoft (<0.5keV) emission 
traces a nuclear outflow cone (Istorchi-Bergmann k, Bonattc 



199UlHielm fc Lindbladlll996l ; lLindbladlll999l ; IVeilleux et al 
2003]). Throughout this paper we assume a luminosity dis- 
tanc<0 of 21 Mpc to NGC 1365. All quoted flux and lumi- 
nosity uncertainties are 90% confidence limits, unless noted 
otherwise. 



Figure 2. Logarithmically spaced X-ray contours and images during the 2002 Chandra observation in the ultrasoft (0.2-0.5 keV; top 
left), the soft (0.5— 1.5keV; bottom left), and hard band (1.5— lOkeV; bottom right). The highest level contours are selected to exclude the 
AGN point sources in each band, while lowest level contours are selected to exclude background emission. The top right panel overplots 
the soft (red), medium (green), and hard band (white) emission contours. The cir cumnuclear star-form ing ring of NGC 1365 is also shown 
here tr aced by the optical super-star clusters (the SSCs, denoted by red squares; Kris ter ; et aL 1993, KJL97 ), the radio hotsplots (green 
circles; IStevens et alJl99Sl SFN99), and the 12 CO molecular hotspots (yellow ellipses; ISakamoto et aTll2007l, S HMMP07). The soft band 
con tours are extende d in the direction of the radio jet (radio hotspot F), and the [O III] A5007 outflow cone of lHielm fc Lindbladl l|l99rj l 
and lLindbladl dl999T ). 



2 THE EXTENDED X-RAY EMISSION OF 
NGC 1365 

2.1 Observation Summary 

NGC 1365 was first imaged by Chandra in 2002, December. 
The observation shows diffuse extended nuclear emission in 

2 NASA/IPAC Extragalactic Database, 

| http : / / nedwww. ipac . ca ltech . edu/ 



the inner 20" x 20" 2 x 2 kpc) an d 23 off-nuclear point 
sourc es in the central ~ 6' x 6' region l|Strateva fc Komossal 
120091 . hereafter SK09). Between 2004 and 2006, NGC 1365 
was the target of extensive X-ray campaigns which resulted 
in six additional Chandra and five XMM-Newton observa- 
tions (with EPIC-pn exposures of 5.4 ks to 56.5 ks); for more 
details of the exposure times and dates of these observations, 
see Table [2] We processed all seven Chandra observations 
with CIAO v.3.4. The 14.6 ks Obs.ID3554 observation is 
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Table 1. NGC 1365: X-ray Observation Summary. T e ff is the 
effective exposure time after corrections for flaring and chip 
position. 



ObsID MJD T off [10 3 s] 




Chandra 3554 


52632 


9.9 (ACIS-S3) 


Chandra 6871 


53835 


13.4 (ACIS-S3) 


Chandra 6872 


53837 


14.6 (ACIS-S3) 


Chandra 6873 


53840 


14.6 (ACIS-S3) 


Chandra 6868 


53843 


14.6 (ACIS-S3) 


Chandra 6869 


53845 


14.6 (ACIS-S3) 


Chandra 6870 


53848 


15.5 (ACIS-S3) 


XMM 0151370101 


52655 


14.4 (EPIC pn) 17.6 (MOS) 


XMM 0151370201 


52679 


2.5 (EPIC pn) 5.4 (MOS) 


XMM 0151370701 


52864 


6.6 (EPIC pn) 8.2 (MOS) 


XMM 0205590301 


53021 


49.0 (EPIC pn) 56.5(MOS) 


XMM 0205590401 


53210 


30.0 (EPIC pn) 44.6 (MOS) 
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Figure 3. Radial profile of the X-ray emission of NGC 1365 from 
2002. The AGN point source is clearly visible in the inner ~ 1 — 2"; 
the X-ray emission associated with the inner spiral extends to at 
least 20". 



reduced to 9.9 ks after the exclusion of a 4.7 ks flare. The 
remaining 6 Chandra observations are free of flaring, with 
effective times ranging between 13.4 and 15.5 ks. 

Using both Chandra and XMM-Newton we can study 
in detail the extended emission of NGC 1365, through both 
spectroscopy and spatially resolved imaging, excluding the 
absorbed and highl y-variable nu c lear po int source described 
in detail by e.g., iRisaliti et all |2005al ). In this paper we 
present visual results based on our original 2002 Chandra 
exposure, while using the full Chandra and XMM-Newton 
datasets in the analysis. 



0.5 1 2 5 

OBSERVED-FRAME ENERGY (keV) 

Figure 4. Chandra spectra of the AGN point source (< 2"; 
rombs), the extended emission (2" < r < 40", see text; trian- 
gles) and the total emission (crosses). The AGN emission dom- 
inates above ~ 2.5 keV, the extended emission dominates below 
~ 1.5 keV. 

2.2 The Inner Spiral and Outflow Cone of 
NGC 1365 

Kinematically, the nuclear region of NGC 1365 shows two 
distinct parts: (1) the rotation field of the inner spiral, and 
(2) the outflowing cone component. The inner spiral is well 
traced by the optical/infrared (IR) continuum, as well as low 
excitation lines like the Balmer lines, and the forbidden [N 

II] , [S II], and [O II] lines. The inner spiral has two promi- 
nent dust lanes (clearly visible on the HST image in Fig. [1]), 
which can be traced b y dense molecular gas, e.g., amonia 
(NH 3 ; lOtt et al.ll2005h . The outflowing cone component is 
traced by higher excitation lines like [O III] AA5007,4958 , 
[Ne III]A3868 and He II A4686 (|Hielm fc Lindbladl Il996l ; 
Lindblad 1999) and can be se e n on Figs. 19 and 22 of 
Lindbladl ifl999l ). IVeilleux et all (|2003h constructed a spa- 
tially resolved map of the narrow-line-region flux ratios ([O 

III] /H/3 and [N II] /Ha) of the central few kiloparsecs of 
NGC 1365. They concluded that starburst-like photoioniza- 
tion dominates the nuclear spiral, with the exception of the 
region associated with the strongest [O III] emission, which, 
together with the region outside 2 kpc, has the characteristic 
ratios of AGN photoionization. 

2.3 Spatially- Resolved X-ray Emission 

The left panel of Figure [T] gives a panoramic view of the 
NGC 1365 spiral, where we have outlined the nuclear region 
that we study in detail. The top right panel of Figure[T]shows 
the 0.3-10 keV X-ray emission contours from the 2002 Chan- 
dra exposure overlaid on the F555W WFPC2 image of the 
inner spirajf]. The 0.3-10 keV X-ray emission traces well the 
nuclear spiral component, suggesting that the extended X- 
ray emission is associated with the circumnuclear starburst 
ring with projected size ~ 20" x 10" (~2 kpex 1 kpc). The 

3 For the multiwavelength comparison we have applied a ~ 1" 
correction to the X-ray positions requiring that the (>0.5keV) 
X-ray emission peak coincides with the position of the optical 
emission peak and the radio center N. 
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Figure 5. Best fit 0.3-10 keV spectrum for Chandra Obs. ID 3554. 
The model fit includes a thermal plasma (apec in XSPEC) with a 
temperature of ~0.7keV and sub-solar abundance, and a power- 
law tail with a photon index T = 1. The two model components 
are shown with dotted lines in the top panel, the fit residuals are 
shown in a separate panel below. 



nuclear starbur st ring; is traced comple t ely by radio 3 cm an d 
6 cm emission (|Forbes fc Norrislll998l ; IStevens et aT] 1 19991). 
and partially by molecular CO emission (S akamoto et al.l 
20071) and optical ly detected super-star clusters (SSCs; 
Kristen et al.lll997l ). as shown in the bottom right panel of 
Figure [T] 

The spatial distribution of the hard (>1.5keV) and ul- 
trasoft (<0.5keV) Chandra X-ray photons provide further 
clues to the origin of the extended X-ray emission. The 2"- 
smoothed X-ray images and contours of the extended emis- 
sion from the 2002 exposure are presented in Figure [5] in 
three different energy bands. The 0.2-0.5 keV X-ray emis- 
sion shown in the upper left panel is not spatially associated 
with the circumnuclear ring; instead it traces the outflow 
cone seen in higher optical excitation lines like [O III]. The 
ultrasoft X-ray emission also coincides with the direction of 
the probable nu clear jet (extending from the nucleus, N, to 
radio hotspot F; IStevens et al.lll999T ). peaking 0.8" SW from 
the centers of the soft- and hard-band emission (which coin- 
cide with the radio center N by construction; top right panel 
of Figure [2} . Following the spatial coincidence between the 
radio, [O III], and ultrasoft X-ray emission, which is con- 
firmed by the emission maps of the remaining 6 Chandra 
observations, we suggest that the ultrasoft emission is ei- 
ther part of an X-ray jet or, more likely, originates in the 
wide-angle outflow. 

The hard (1.5-10 keV) band emission shown in the lower 
right panel of Figure [2] coincides with some of the radio and 
molecular hotspots and SSCs. The converse is not true - 
not all the SSCs have hard X-ray counterparts - the bright- 
est SSCs (S10 and S3), for example, do not correspond to 
hard X-ray peaks. The general impression on this and the 
remaining 6 Chandra exposures of the same region is that 
the hard X-rays trace well the circumnuclear starburst ring 
vissible in the optical, radio, and sub-mm. This further sug- 
gests that the extended hard X-ray emission is associated 
with the vigorous nuclear star formation. 

The soft (0.5-1.5 keV) X-ray emission is shown in the 



bottom left panel of Figure [2] It traces the inner spiral, em- 
phasizing the northern spiral arm in contrast to the adjacent 
dust lane. The radio hotspots D, E, and G as well as their 
12 CO molecular counterparts fall within this dust lane and 
are absorbed in the soft X-rays. The southern spiral arm 
is much fainter and not visible in the X-ray band, because 
the NW part of the spiral is nearer to us, while the SE arm 
is seen absorbed through the galactic plane (see also Fig- 
ure [TJ. T he soft and hard X-ray emi ssion peaks coincide in 
Figure l2l Irtisaliti et all |2005al , f2007l ') showed that the hard 
X-ray emission comes from the central AGN, which is always 
obscured by a neutral column of at least ~ 2 x 10 23 cm -2 . 
This means that no direct soft X-ray emission is expected 
from the AGN below 1.5 keV, implying either indirect AGN 
emission (scattering) or an extra unresolved starburst com- 
ponent which we discuss below. Considering the observed 
PSF profile shown in Figure [3] the contributions from the 
wings of the PSF cannot dominate the observed soft emis- 
sion from the center. 



2.4 Chandra spectra of the extended emission 

Chandra's spatial resolution allows us to model the ex- 
tended emission of NGC 1365 by excluding the AGN point 
source. As can be seen on Figure [3] which shows the radial 
profile of NGC 1365 during Obs. ID 3554, the AGN point 
source dominates the emission from the inner ~ 2", which 
we excluded from the following analysis. Figure U shows 
the Chandra 2002 spectra extracted from the inner 2", a 
40" circular annulus excluding the inner 2" (the extended 
emission), and the total emission within 40". The emission 
above ~ 3keV is clearly dominated by the active nucleus, 
that below ~ 1.5 keV - by the exte nded emission, as al- 
ready noted by iRisaliti et al.l (|2005aT ). The extended emis- 
sion peaks around 0.7-0.8 keV, which suggest a thermal- 
plasma emission origin. 

After experimenting with a few different thermal- 
plasma models (XSPEC models mekal, cmekal, apec, ray- 
mond, etc.), we find that apec models the thermal emis- 
sion best, but the differences are not significant^ In the full 
0.3-10 keV Chandra band, in addition to the thermal-plasma 
component, there is an extended power-law emission tail 
above 2 keV. Using the 2-10 keV extended-emission spec- 
trum, we determined a best-fit photon index ofP = 1.1±0.5 
in ObsID 3554, which is substantially different from that of 
the AGN in the 2-10 keV band, whose simple power law 
slope is F w —0.8; this suggests that the extra compo- 
nent is not dominated by emission from the AGN0 The 
addition of the power-law component significantly imroves 
the 0.3-10 keV fits, with F-test probability of 10" 8 . The 
best fit 0.3-10 keV spectral model including the thermal 



4 All models include absorption equal to the Galactic value of 
1.4 X 10 20 cm" 2 . 

5 In our original 2002 observation, the extended emission counts 
in the 2-10 keV (5-10 keV) band are equal to about 28% (56%) 
of the total AGN counts in the same band, while the expected 
contamination from the wings of the PSF should be about 10%. In 
the subsequent Chandra observations, when the AGN was more 
luminous, the contamination is stronger, especially around and 
above 5 keV, which is also confirmed by comparing the hard-band 
spectral slope of the extended emission with that of the AGN. 
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Table 2. NGC 1365 0.3—1.5 keV emission: Fluxes and Spectroscopic Fits. All model fits include a cold absorber equal to the Galactic 
value of 1.4 X 10 20 cm -2 ; the second fit to ObsID 3554 includes intrinsic absorption of (7± 2) X 10 20 cm -2 . The fluxes are quoted in 
units of 10 — 13 ergs -1 cm -2 and all XMM-Newton fluxes are scaled by 20% to account for the AGN contribution. All errors are 90% 
confidence limits; fit parameters given in parenthesis were held fixed during the fit. Approximate values are given for parameters 
which were unconstrained by the fit. 



ObsID 


Counts 


model 


-^0.3-1.5 koV 


r 


fcTj [keV] 


Zi/Z B 


kT 2 [keV] 




x 2 /DoF 


Q tr C /I 
O004 


1 ^oi 
loy i 


apec-{- zpow 


c q+0.6 




fi fix. _|_ n HQ 
U.DO it U.Uo 


ni7J.fl C\A 
U. 1 ( it U.U4 






4o/ 04 


Do I 1 


ZUUD 


apec-\- zpow 


5 - 4 -i.l 


( L ) 


n K/i -U n no 
U.04 it U.UZ 


n 1 « J- n nQ 
U. 10 it U.Uo 






A Q i A fi 
4o/4D 


Do ( Z 


A A I 


apec-\- zpow 


5 ' 5 -l .4 


( L ) 


n fi/i -U n no 

U.04 it U.UZ 


n oi J- n no 

U.Zl it U.UZ 






f^n /f^fi 

OU/ OD 


Do f o 


Ol QO 


apec-\- zpow 


5 ' 2 -1.4 


K L ) 


U.DO ± U.Uo 


U. 10 it U.U4 






If / 00 


DoDo 


OOOI 
ZZZ1 


apec-\- zpow 


5 ' 4 -1.0 


I 1 ) 


n 71 j. n nQ 
U. f 1 it U.Uo 


n i k J- n nQ 
U. 10 it U.Uo 






Dy/ oy 


DoDy 


ZoOO 


apec-\- zpo w 


tr 4+0.5 
5 ' 4 -1.3 




n -U n no 
U.Do it U.UZ 


n on J- n n/i 

U.ZU it U.U4 






QQ /c;^ 
OOf Do 


Do / U 


zlyo 


apec-\- zpow 


c: c+0.5 
°-°-0.9 




n fio J_ n no 
U.DZ ± U.UZ 


n i k J_ n nQ 
U. lo it U.Uo 






DO/ Oo 


m Q7m m 

UlOlo ( U1U1 


yio ( 


apec-\- zpow 


fi ^+ - 3 

b ' 5 -0.3 


( L ) 


n fi^ _i_ n m 
U.DO it U.U1 


n i n J- n m 
U. 1U it U.U1 






Q/1 fi/QIO 
o4o/olo 


m f^i Q7nom 

UlOlo < UZU1 


OflQfl 
ZUoU 


apec-\- zpow 


fi fi+ ' 8 


( L ) 


n vn -U n nQ 
U. (U it U.Uo 


n no _i_ n m 
u.uy it U.Ul 






OA /Q7 

yo/o * 


U -LU -LO 1 U 1 \J _L 


4575 


{t'rifif 1 ynntn 


6 8+ ' 4 


(i \ 

X 1 ) 


n 7i + o m 

\J • 1 -L _1_ U.Ul 








191/175 


0205590301 


35873 


apec+zpow 




(1) 


0.64 ±0.01 


0.096 ± 0.004 






702/626 


0205590401 


22517 


apec+zpow 




1.0 ±0.5 


0.64 ±0.01 


0.10 ±0.01 






542/468 


3554 


1591 


zwabs *apec 


7 7 +0.0 
'■'-1.3 




0.66 ±0.03 


0.07 ±0.01 






52/54 


6871 


2006 


apec+ apec 


5 5+ - 4 
°'°-0.8 




1.2 ±0.2 


0.06 ±0.03 


0.61 ±0.04 


(1) 


38/45 


6872 


2276 


apec+ apec 


5 4+ ' 4 




1.0±0.1 


(0.1) 


0.55 ±0.04 


(1) 


60/56 


6873 


2182 


apec+ apec 


5 4+ ' 4 

-0.7 




0.99 ±0.07 


0.10 ±0.02 


0.54 ±0.07 


(1) 


63/54 


6868 


2221 


apec+ apec 


^ 4+0.6 




1.7 ±0.4 


0.70 ± 0.03 


0.70 ± 0.03 


(1) 


68/58 


6869 


2366 


apec+ apec 


5 4+ ' 5 




1.2 ±0.2 


0.10 ±0.03 


0.58 ±0.04 


(1) 


73/54 


6870 


2193 


apec+ apec 


5 4+ ' 7 

J * -0.7 




1.2 ±0.2 


<o.i 


0.62 ±0.03 


(1) 


68/52 


0151370101 


9167 


apec+ apec 


6 6 +0 ' 4 

D.D_ 2 




1.5 ±0.1 


<0.02 


0.66 ±0.01 


(1) 


350/317 


0151370201 


2080 


apec+ apec 


6 7+ ' 6 
"•'-0.7 




1.3 ±0.2 


< 0.02 


0.66 ±0.03 


(1) 


99/86 


0151370701 


4575 


apec+ apec 


6 7+ ' 4 

"•'-0.5 




1.1 ±0.1 


0.02 ±0.01 


0.66 ±0.02 


(1) 


190/174 


0205590301 


35873 


apec+ apec 


7 q+0.2 

'•"-o.i 




1.7 ±0.1 


< 0.02 


0.64 ±0.01 


(1) 


724/625 


0205590401 


22517 


apec+ apec 


6 9 +0 1 
"• M -0.2 




1.15 ±0.05 


0.020 ± 0.005 


0.62 ±0.01 


(1) 


547/468 



plasma and power-law emission components for Chandra 
ObsID 3554 is shown in Figure [5] The 0.3-10 keV luminosi- 
ties of each of the two components are: 3.0^0 2 x 10 40 ergs -1 
for the thermal component and l.Si^ii x 10 40 ergs -1 for the 
power-law component. As noted above, the thermal plasma 
dominates the soft-band extended emission; in the 2-10 keV 
band the power law model contributes 1.4jlj"| x 10 40 ergs -1 
to the extended emission, while the thermal plasma only 
~ 1.3 x 10 39 ergs -1 . 

For some of the remaining Chandra observations, e.g. 
ObsID 6868, the 2-10 keV emission seems to be contami- 
nated by the AGN even beyond 2" from the center, as evi- 
denced by the 2-10 keV power law photon index of V ~ —0.5, 
the large ratio between AGN and extended-emission counts 
around 5 keV, and the residuals between 6 and 7 keV (which 
are consistent with an Fe-Ka origin). Consequently, we re- 
strict our analysis to the 0.3-1.5 keV band, which, as shown 
above both spatially and spectroscopically, is dominated by 
the extended emission, even for the case of the much lower 
spatial resolution XMM-Newton data. Table 12.41 lists the 
total counts and thermal plasma model fit parameters ob- 
tained in the soft band for all seven Chandra observations. 
The single temperature thermal plasma model provides a 
poor fit to the 0.3-1.5 keV data in all cases, and we consider 
the addition of a power law or a two-temperature plasma 
model instead. Both models provide significant improve- 
ments of the spectral fits (F-test probabilities of < 0.1%) 
and in most cases the current data does not show a clear 



preference for either of the two models. For the thermal 
plasma plus power-law model results shown in the top part 
of Table [2~4l we fix the power law slope to F = 1, the value 
found in the 2-10 keV Chandra fits unaffected by the AGN. 
This slope is confirmed in the soft band fit to the longest 
XMM-Newton observation, as discussed below. An exam- 
ple Chandra 0.3-1.5 keV spectral fit from the upper part of 
Table EOKObs. ID 6870) is shown in Figured 

Overall the Chandra spectra of NGC 1365 suggest a 
thermal plasma extended emission origin for the soft X-ray 
emission, with additional contribution from a second tem- 
perature thermal plasma or a F = 1 power law. The fit- 
ted abundances in Table [2~4l are sub-solar, inconsistent with 
our general ideas about the recycling of materials around 
AGNs and the abundance measurements and gradients of 
NGC 1365 (e.g., Galliano et al. 2005 report central abun- 
dance of 3 times solar). Sub-solar abundances are commonly 
found when applying simple thermal plasma models (e.g., 
iBuote fc Fabianl Il998l ). and are considered a sign of the 
model simplifications (e.g. thermal equilibrium vs. shocks, 
etc.). 

From Figure [4] the spectrum of the inner 2" emis- 
sion has a low flux soft X-ray component extending be- 
low 2.5 keV. As noted in § 12.31 this unresolved core emis- 
sion in the soft band is unlikely to originate directly from 
the active nucleus, which is absorbed by a large neutral 
column. The 0.3-1.5 keV core (within 2") emission can be 
represented by a spectral model including thermal plasma 
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Figure 6. Best fit 0.3-1.5 keV spectra for an example Chandra 
(Obs. ID 6870, top) and XMM-Newton (Obs. ID 0205590401, bot- 
tom) observation. The model fits include a thermal plasma (apec 
in XSPEC) with a temperature of ~0.62keV, and a power-law 
tail with a photon index T = 1. The model fit components are 
shown separately with dotted lines below the fit; the fit residuals 
are shown in the lower panel. 



(kT = 0.8 ± 0.1 keV and abundances fixed to 0.17 of solar) 
and a power law absorbed by a Galactic neutral absorber 
consistent with the model fit of the extended emission. The 
luminosity of this core component is low - L . 3-1.5 keV = 
4.4±1.1 x 10 39 ergs s _1 , w 15±5% of the 0.3-1.5 keV band lu- 
minosity of the extended emission, too large to be attributed 
to the PSF wings of the extended emission surrounding the 
center. Based on the similarity of the spectral shape of this 
core component to that of the extended star-formation asso- 
ciated emission, we speculate that this soft X-ray core com- 
ponent is dominated by unresolved recent star formation in 
the immediate vicinity of the active nucleus. 



2.5 XMM-Newton spectra of the extended 
emission 

Due to the fact that the soft emission is dominated by 
the extended component of NGC 1365 and the heavily ob- 
scured AGN's contribution is expected to be negligible be- 
low 1.5 keV, we can use the 0.3-1.5 keV XMM-Newton ob- 
servations of NGC 1365 to constrain the properties of the 
extended emission. The unresolved soft X-ray emission com- 
ponent found within 2" of the active nucleus in the Chan- 
dra observation above, has identical 0.3-1.5 keV spectrum 
and will only increase the total flux by about 15%. In 
each case we extract the 0.3-1.5 keV XMM-Newton spec- 
trum from a 40" region centered on NGC 1365 and bin it to 
>20 counts per bin in order to use \ 2 minimization fitting 



in XSPEC. The corresponding XMM-Newton spectral fits 
are given in Table 12.41 together with the Chandra fits for 
comparison. The quoted fluxes are XMM-Newton pn fluxes. 
Note that they are about 20% brighter than the correspond- 
ing Chandra fluxes, in agreement with the expected soft- 
band contribution of the unresolved core emission. The 0.3- 
1.5 keV XMM-Newton spectra are well represented by the 
thermal plasma plus power-law model found in the Chan- 
dra fits, with the exeption of the two longest observations, 
ObsID 0205590301 and ObsID 0205590401, where this fit is 
inadequate (x 2 / DoF=702/625 and \ 2 7 '-DoF=542/468, re- 
spectively). Adopting varying abundances for some elements 
(e.g. Fe, Mg, Si, He) provides acceptable fits in both cases 
(x 2 /£>oF=663/621 and x 2 /£>oF=514/464), but consider- 
ing the simplicity of the one-temperature collisionally ion- 
ized diffuse gas models, the caveats associated with the apec 
and mekal models in XSPEC, and the inconsistencies and 
poor constraints on the fitted abundance values of these 
model fits, we refrain from inter preting these abundance 
variations here. IWang et all (|2009l ) suggest that allowing ele- 
ment abundances to vary independently for different spatial 
regions can resolve the problem of the small fitted overall 
abundance. 

We checked the longest XMM-Newton observation for 
varibility in the soft band (0.3-1.5 keV) on timescales of 
hours, but found no significant variation. 



3 DISCUSSION AND CONCLUSIONS 

3.1 Multiwalength View of the Nuclear Spiral 

The spatially-resolved 0.5-10 keV emission of NGC 1365 
traces well the nuclear spiral component, suggesting that 
the extended X-ray emission is associated with the circum- 
nuclear starburst. This is confirmed by the spatially resolved 
Chandra spectroscopy, which, together with the unresolved 
XMM-Newton soft-band spectra, are dominated by thermal 
plasma emission with temperature ~0.7keV. The ultrasoft 
X-ray emission (<0.5 ke V), on the other hand, is c onfined to 
the outflow cone seen bylHielm fc Lindbladl (fl99^ V iLindbladl 
\ 19991 ). and lVeilleux" et al. I (|2003l ) in the kinematic signature 
of the high excitation optical emission lines. Soft-X-ray emis- 
sion associated with a NLR outflow has been seen before, no- 
tably in NGC 4151 and NGC 1068, where the o utflows are 
photoionized and photoexcited b y the AGN fe.g-. lQgle et alj 
200d. 120031; ISchurch et ai1l2004l ). In the case of NGC 1365, 
Veilleux et al.l l|2003l ) show that the optical narrow-line ra- 
tios associated with the outflow cone are consistent with 
AGN photoionization outside the nuclear spiral region, but 
have the photoionization characteristics of a starburst inside 
the inner 1-2 kpc. In fact, as we argue below, the NGC 1365 
starburst is powerful enough to drive a strong outflow like 
the one observed. 

The molecular starburst ring of NGC 1365 is de- 
fined in the radio band, by radio 3 cm an d 6 cm emission 
(|Forbes fc Norrislll998l ; IStevens et al"1ll999T I and well traced 
in the MIR by Spitzer, as well as by the string of su per-star 
clusters in the optical (SSCs; iKristen et alj [l997h . There 
is evidence for molecular gas in the center which is dense 



(> 10 cm ) and cold ( kinetic temperature of 40 K for 
NGC 1365; lOtt et al.ll2005l ). traced by NH 3 , which correlates 
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Figure 7. Spitzer IRAC3 (5.8 /an, left) and IRAC4 (8.0 (im, right) images of the nuclear starburst ring. The central nuclear emission 
dominates at and below ~6 /im; the circumnuclcar starburst ring emission is comparable at 8 (im, and dominates the IR emission at longer 
wavelengths, for example, the 24 /im MIPS emission shown in the left panel as contours. The right panel displays also the 0.3-10 keV 
Chandra contours (from Obs. ID 3554). The super-star clusters, radio hotspots, and 12 CO molecular hotspots of the nuclear starburst 
ring are shown in the right panel with symbols following those of Figure [T] 



well wit h other star-formation indicators. The CO emission 
maps of ISakamoto et all (|2007l ) confirm the presence of less 
dense molecular gas, with the 12 CO and 13 CO molecular 
maps tracing the circumnuclear ring structure dotted by 
emission hotspots. 

NGC 1365 wa s recently imaged w ith Spitzer with both 
IRAC and MIPS (|Brandl et alj|2006h . Figure [7J shows the 
5.8 /im and 8 /im IRAC images together with the 24 /im 
MIPS (left panel) and the 0.3-10 keV X-ray (right) contours, 
as well as the circumnuclear ring SSCs, radio, and molecular 
CO hotspots. The circumnuclear star forming ring is trace- 
able in the MIR in both panels, gaining prominence with 
respect to the nuclear component at 8 /im. At shorter MIR 
wavelengths the nuclear point source dominates the emis- 
sion, at wavelengths increasingly longer than 8 /im the star- 
burst dominates. The 24 /im MIPS emission contours shown 
in Figure [7] peak in the NW-running dust lane, where the 
X-ray emission is mostly absorbed in the Chandra band, but 
the molecular 12 CO emission (COl, C02, and C03 in Fig- 
ure [JJ and the radio emission (hotspots D, E, and G) are 
strong. We address the to the nature of the emission in this 
strong MIR emission region in the next section. 

The (spatially unresolved) 5-38 /im spectrum of 
NGC 1365 is a fairly shallow power law with no evidence 
of silicate absorption and weak 6.2, 7.7, 8.6, and 11.3 /im 
polycyclic aromatic hydrocarbon (PAH) features. The star- 
burst dominates, but a prominent AGN component must 
also be present, and overall the nucleus lies (fully on the star- 
burst side but) close to the line sepa rating starbur s ts an d 
AGN on the MIR-slope diagnostic of iBrandl et all {2006). 
Using our assumed distance of 21 Mpc and the Spitzer the 



15 /xm and 30 /im Spitzer fluxes (jBrandl et ail [2006). the 
MIR SFR estimate is SFR Sp ii Ze r ~ 9M© yr" 1 . This is about 
25% smaller than the FIR-estimated global NGC 1365 SFR, 
which is SFRpm ~ 12M Q yr _1 , if we assume the Lonsdale 
et ai 1985 Lfir = 6.8 x 10 10 Z/© measurement and eqn. 3 of 
iKennicuttJ (| 19981 ). 

The soft-band X-ray luminosity can also be used to esti- 
mate of the star-for mation rate. With the help of eqn. 14 of 
iRanalli et ail (|2003l ). SFR = 2.2 x 10 -40 Lo.5-2kcv M© yr -1 , 
and the observed 0.5-2 keV luminosity of the extended emis- 
sion, L .5-2kcv = 3 — 3.6 x 10 40 ergss _1 (where the higher 
value includes the contribution of an unresolved starburst 
component within 2" from the AGN), we estimate a soft-X- 
ray of SFR .5_2kov = 7 — 8M Q yr~ . This is in good agree- 
ment (10-20% lower) than the Spitzer nuclear SFR estimate, 
which covers a slightly larger area, and about 30-40% lower 
than the global FIR-estimated SFR. We conclude that most 
of the strong star formation (60-70% or 7-8 Mq yr -1 ) is con- 
strained to the nuclear starburst ring in the inner region of 
NGC 1365. 

3.2 Radio hotspots, Super-Star Clusters and 
Molecular clouds 

Radio hotspots A, D, G, and H have flat, non- t herma l 
spectra betwee n 6 a nd 20 cm (|Sandqvist et ail Il995l ). 
ISandqvist et ail (|l995l ) suggest that radio hotspots A, D, 
and G (A and G are unresolved at 2 cm, while D is only 
marginally detected) are related to the star formation and 
supernova activity. Radio hotspots D and G (which co- 
incide with CO hotspots COl and C03) fall within the 
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Figure 8. X-ray spectra of select multi-wavelength counter- 
parts for Chandra Obs. ID 3554: the radio-hotspot A (SSC10; top 
panel), radio hotspots D&G (middle panel), and radio hotspot 
F (jet/high- ionization outflow; bottom panel). 



prominent dark NW dust lane, and are associated with 
tw o obscured star- forming regions detected in the infrared 
by iMoorwood et al and visible in the IRAC image 

of N GC1365 shown in iBrandl et all (120061) and our Fig- 
ure [3 [Galiiano^eral] (|2005h and lSakamoto et alj (|2007T ) sug- 
gest that the MIR peaks associated with radio hotspots 
D, E , and G (M4, M5, and M6 in iGalliano et "ail 120051 . 
120081) correspond t o SSC s. Using VLT images and spec- 
tra IGalliano et all (|2008h revised the mass and age esti- 
mates of earlier works to show that the superstar clusters 
in NGC 1365 are the most massive clusters found so far, 
with masses of at least 10 7 Mq and ages of ~ 7Myr. These 
SSCs are still embedded in large amounts of gas and dust 
due to dust trapping, unlike typical lower-mass star clusters 
found in less dense environments which sweep away their 
gaseous/dusty envelopes much earlier. 

Using the X-ray images we can obtain estimates of the 



X-ray radiation from the radio and 12 CO hotspots and the 
HST SSCs. X-ray spectra extracted from l"-region^f| around 
the radio positions of sources D and G have only a handful 
of counts (16-24), but are consistent with the apec+zpow 
(or simple apec) models with parameters fixed to those re- 
ported in Table 12.41 for each Chandra observation, if we 
also include an intrinsic absorber with column density of 
~ 0.5 - 1 X 10 22 cm" 2 (see the middle panel of Fig. |5J. The 
unabsorbed luminosity estimates for the 0.3-10 keV band 
correspond to L x ~ 1 — 3 x 10 39 ergs -1 for both radio 
hotspots taken together. These luminosities are in agree- 
ment with the theoretical p redictions fo r the highest mass 
clusters (10 6 M©) studied bv lSilich et"al] (120051 ). which have 
L x ~ 1 — 10 x 10 39 ergs -1 (depending on the density of the 
surrounding interstellar medium, the gas metallicity, etc.), 
but these predictions do not take into account the effects of 
gas-trapping in heavier clusters, which is likely relevant for 
the SSCs behind radio hotspots D and G. 

Radio hotspot A, which has no molecular 12 CO equiv- 
alent, falls in an optically and X-ray unabsorbed region and 
is likel y associated with a s uper-star cluster SS C10 identi- 
fied bv lKristen et all l|l997l ) (|Stevens et al.lll999T l. There are 
15-39 counts within 1" of the optical (similar for the slightly 
offset radio) positions in the different Chandra observations 
and apec+zpow Cash statistic model fits with parameters 
fixed to those of the extended Chandra emission fits from 
Table 12.41 and free normalization provide acceptable fits in 
all cases (see the top panel of Fig. [5] for an example). The 
inferred 0.3-10 keV 1" flux is £ x ~ 6 - 10 x 10 38 ergs _1 . 
The B-band and radio luminosities of SSC10 are L B = 4.2 x 
10 39 ergs -1 and L 3cm = 5.2 x 10 36 ergs -1 , the corresponding 
ratios L X /L B ~ 0.1 - 0.2 and L 3cta /L* ~ 0.005 - 0.009. 

Radio hotspot F (which coincid es with CQ5) is tradi- 
tionally ident i fied w ith a radio jet l|Sandqvist et al.l 1 19951 : 
IStevens et all I1999T I. We argued above, that the ultra- 
soft X-ray emission is associated with the outflow cone 
seen in higher excitation lines like [O III] and [Ne III] 
l|Hielm fc Lindbladlll996l : lLindbladlll999l) . We extracted X- 
ray spectra from the region between radio hotspots N and F, 
which overlaps with the ultrasoft emission contours shown 
in Figure [2] for each of the Chandra pointings, excluding 
the 2"-radius around the center which is affected by the 
AGN. The spectra can be fit by apec+zpow models, with the 
apec model parameters fixed to these quoted in Table 12.41 
as shown in the bottom panel of Fig. [8] Unlike the case 
of the overall extended emission component however, the 
power-law component dominates the X-ray spectrum of the 
outflow at both high and low energies, i.e. the power-law 
index is much steeper - typically V ~ 2.5 - and the nor- 
malization of this component much higher. The apec+zpow 
fits (including absorption equal to the Galactic value in all 
cases) are acceptable in most cases, but the residuals show 
excess emission below 0.5 keV, i.e. the ultrasoft emission 
seen in Figure [2] is not fully accounted for by the steeper 
power-law. The lower limit to the energy (since we exclude 
the highest surface brightness inner region due to AGN 
contamination) in the 0.3-10 keV band for this outflow is 
Lxjet > 4 x 10 39 ergs -1 . 



e For Chandra on-axis observations, 80—90% of the total energy 
falls in a circle with radius ~ 1". 
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3.3 Starburst-Driven Superwind 

Starburst-driven superwinds (as well as AGN driven out- 
flows) can explain the extended, often biconical. X-ray 
emission in nearby gal ax ies (e.g.. [Heckman et all 1 19961: 
iKomossa fc Schulzl 1 19981 ). iKomossa fc Schulzl (|l99a ) con- 
cluded that the extended soft X-ray emission observed in 
NGC 1365 co uld be the result of a sta rburst driven wind. 
We repeat the IKomossa fc Sch ulz ( 1998) estimates here, us- 
ing the new data on the SFR and new supernova (SN) rate 
estimates, assuming a velocity of ~100-200 kms -1 on ~1— 
2 kpc scales. 

We use two different methods to obta in an estimate of 
the SN rate. ICappellaro fc Turattol (|200ll ) tabulate the ex- 
pected SN rates as a function of B-band luminosity and 
the galaxy type. The B-ba nd luminosity of NGC 1365 is 
7.9 x 10 10 L© (see Table 2 of ISandqvist et all fl988l ). result- 
ing in a SN rate of 0.09 yr _1 . Using the FIR luminosity 
method, w e obtain slightly higher SN rate in the range 0.14- 
0.16 yr~ llvan Buren fc Greenhouse! [l99l; |m 

annucci et al.l 

l2003ft . Both the B-band and the FIR-luminosity SN rate 
estimates apply to the whole galaxy. Scaling the global SN 
rate by assuming that 60-70% of the new stars and SN are 
born in the inner spiral, we arrive at a central SN rate in 
the range 0.05-0.1 yr _1 . For comparison the radi o-band in- 
ferred SN rate reported by IStevens et all l|l999h is about 

0. 02-0.03 yr~ x over a smaller (~lkpc 2 ) area including the 
central ~0.6kpc. 

Following IKomossa fc Schulzl ( 19981 ). the mechanical 
power or a starburst driven wind in NGC 1365 is L moch ~ 
1 — 3 x 10 42 ergs -1 , assuming each supernova gives L SN = 
10 51 ergyr _1 . With a soft-band luminosity of L0.1-2.4kcv 

1, ; x 10 40 ergs _1 , the IKomossa fc Schulj ll 19981) estim ates 
(their eqns. 1 through 3, based on lHeckman et a.1.11 1996T ) still 
hold, if we assume a mean density of the swept-up interstel- 
lar medium in the range n — 0.1 — 1cm -3 for an outflow 
velocity of ^100-200 kms -1 on -~l-2kpc scales. Therefore 
the observed starburst could supply the observed thermal 
X-rays in a wind driven shell. 

The Ha emi ssion expected fo r a SFR of ~7-8M0yr _1 
is ~ 10 42 ergs -1 jKennicuttlll998l) much lar ger than the Ho 
luminosity observed bv lKristen et al.l (| 1997T ) : 8 x 10 40 ergs -1 
(converted to a distance of 21 Mpc) for the two most lumi- 
nous Ha hotspots located SW of the nucleus (in the region 
of the HST-detected SSCs shown in Figures [T] and 0. 
The total extended Ha luminosity likely exceeds this esti- 
mate by a factor of a few. Considering the large obscuration 
inferred towards the star-formation regions in the NE of the 
nucleus (for example, those associated with radio hotspots 
D, and G, which are optically faint, but infrared and radio 
bright, hidden behind Nh ~ 10 22 cm -2 according to the X- 
ray emission), it is conceivable that as little as ~10% of the 
Ha luminosity associated with the vigorous star formation 
is observed directly. 



3.4 AGN contribution and IR-Optical- X-ray 
relations 

IKomossa fc Schulzl (j 1998h considered the place of the 
NGC 1365 AGN among t h e IR -Optical-X-ray rela- 
tions found by I Ward et al.l (|l988l ) for the sample of 
iPiccinotti et al.l (|l982l ). They concluded the AGN was too 



faint in the X-rays relative to the luminosity expected 
based on the observed broad Ha-line emission and the 
total IR emission. We know now that t he AGN lumi- 
nosity in the 2-10 keV band, based on the llvomoto et al.l 
(|l997h fit to the ASCA data, was underestimated, and 
that in the Compton thin state, the AGN absorption- 
corrected luminosity is about an order of mag nitude higher, 
L AGN 2-10 key = 0.8 - 1.7 x 10 42 ergs _1 l|Risaliti et all 
l2005al . 12007ft . The value expected for the hard-band X-ray 
lu minosity base d on th e broad-line Ha luminosity observed 
by ISchulz et~aH (| 19941 ) . is L A GN,2-iok C v ~ 0.6 x 10 42 ergs _1 
(corrected to a distance of 21 Mpc), comparable to the 
lower range of the observed hard-band AGN emission. 

Usin g the IR con ti nuum fluxes at 15 /im and 30 )im 
given by iBrandl et al.l (|2006l . see their Table 3), we esti- 
mate L 25 - 6 o^m = 2.3 x 10 44 ergs _1 (as in I Ward et ailll988l . 
including setting Ho = 50 kms -1 Mpc -1 for the compari- 
son). The expected hard-band X-ray lumino s ities is then a 
factor of 10-20 higher than the lRisaliti et all l|2005bl ) obser- 
vations. We believe this discrepancy is primarily due to the 
fact that the starb urst domina t es th e IR emission beyond 
6-15 fj,m, while the lWard et alj l|l988l ) relations refer to the 
ha rd-band X-ray and IR emission of the AGN. According 
to I Ward et~ai1 (|l988l ). the IR emission from the AGN com- 
ponent is NGC 1365 should be L^-eo^m < 2 x 10 43 ergs _1 . 
If the AGN dominates the emission at and below 6 /im, the 
6 fim flux density should be a better predictor of the AGN 
hard-band flux. This is indeed the case. The average ra- 
tio of the 6 /jm flux density to the 2-10 keV flux for type 2 
AGNs was found to be log[F , 2 -i kcv/^^(6/im)] = —0.61 by 
iLutz et al.ll|2004r ). in agreement with the values observed for 
NGC 1365: -0.9 «C \og[F 2 - W ± BV /uF„(6 ftm)] < -0.6. 



4 SUMMARY 

Using Chandra's excellent spatial resolution and XMM- 
Newton's higher photon gathering capabilities we study the 
extended nuclear emission of the supergiant barred galaxy 
NGC 1365. The combination of radio, infrared, optical and 
X-ray data reveals the complex structure of the nuclear re- 
gion (within about a kiloparsec of the center) highlighting 
the competing effects of the AGN, the circumnuclear star- 
burst ring, and the central outflow. 

The ultrasoft X-ray emission (below ~0.5keV) 
is spatially coincident with the AGN-driven outflow 
cone traced by higher excitation optical emission lines 
ilHielm fc Lindbladl 1 19961 : iLindbladl 1l999l : IVeilleux et ail 
120031) . The X-ray spectrum of the outflow region can be 
represented by a dominant steep (F «2-3) power-law com- 
ponent; the steepness of the power law is driven by the excess 
emission below ~0.5keV. 

Spectroscopic as well as spatial evidence suggests that 
the soft X-ray emission (0.5-1.5 keV) of NGC 1365 has ther- 
mal origin related to the ~1 kpc-radius circumnuclear star- 
burst. Thermal plasma emission models can fit the Chandra 
and XMM-Newton spectra well, with kT ~0.6-0.7 keV; a 
second thermal component with kT ~1.3keV or a power- 
law with r ~ 1 is also necessary to fit the observed spec- 
tra. Considering the high star-formation rate in the nuclear 
region (~7-8 M Q yr _1 ) and the associated high SN rates, a 
starburst-driven superwind is likely present, and can account 
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for the wide angle outflow observed in the higher excitation 
optical emission lines and the ultrasoft X-rays. 

The hard X-ray emission beyond ~2 keV is dominated 
by the obscured AGN, which also contributes substantially 
to the infrared emission below ~6-15 fim. The circumnuclear 
ring is also traced by low-surface brightness hard X-ray emis- 
sion, suggesting that an obscuring column of « 10 22 cm -2 is 
hiding the ~ 7 Myr old (but still shrouded in gas and dust) 
~ 10 7 Mq super-star clusters seen in the radio, molecular 
CO, and mid-infrared in the NE part of the nuclear star- 
burst ring from view in the optical and soft X-ray bands. 

The central regions of NGC 1365 are influenced by the 
strong bar and nuclear spiral, which channel sufficient ma- 
terial to the innermost regions and are likely responsible 
for the circumnuclear starburst as well as fueling the super- 
massive black hole. Overall, the extended emission (which 
also contains unresolved point sources) contributes about 
5 x 10 40 erg s _1 in the Chandra X-ray b and, similar to the to- 
tal em ission from X-ray point sources (|Strateva fe Komossal 
120091 ) during the 2002 Chandra observation. The AGN is 
an order of magnitude more powerful (1 — 2 x 10 42 ergs" 1 ) 
than the combined extended plus point source emission 
(~ 10 41 ergs" 1 ), but ~70% of this emission is obscured by 
gas in the immediate vicinity of the active nucleus, which 
is associated with the broad-line clouds. The moderate- 
luminosity AGN is often outshined by the nuclear starburst 
(e.g., in the radio and 24 /im infrared, as well as soft 0.5- 
1.5 keV X-ray emission), and both the infrared and X-ray 
emission result from the complex interplay of the stellar and 
AGN emission modified by intrinsic absorption, as well as 
absorption related to orientation of the inner spiral with re- 
spect to our line of site. 



ACKNOWLEDGEMENTS 

This project has been partly funded by the DFG Prior- 
ity Programme 1177 "Galaxy Evolution". We thank Peter 
Predehl for useful discussions and a critical reading of the 
manuscript. 

This work is based on observations obtained with 
XMM-Newton, an ESA science mission with instruments 
and contributions directly funded by ESA Member States 
and the US (NASA). In Germany, the XMM-Newton project 
is supported by the Bundesministerium fur Wirtschaft und 
Technologie/Deutsches Zentrum fur Luft- und Raumfahrt 
(BMWI/DLR, FKZ 50 OX 0001) and the Max-Planck So- 
ciety. It is also based on an observation obtained with the 
Chandra X-ray telescope (a NASA mission). 



REFERENCES 

Alloin, D., Edmunds, M. G., Lindblad, P. O., & Pagel, 

B. E. J. 1981, A&A, 101, 377 
Brandl, B. R., et al. 2006, ApJ, 653, 1129 
Buote, D. A., & Fabian, A. C. 1998, MNRAS, 296, 977 
Cappellaro, E., & Turatto, M. 2001, The Influence of Bi- 
naries on Stellar Population Studies, 264, 199 
Chen, P. S., Fu, H. W., & Gao, Y. F. 2005, New Astronomy, 
11, 27 



Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, Q. F., 
Perley, R. A., Taylor, G. B., & Broderick, J. J. 1998, AJ, 
115, 1693 

Donato, D., et al. 2001, ApJ, 375, 739 
Edmunds, M. G., & Pagel, B. E. J. 1982, MNRAS, 198, 
1089 

Edmunds, M. G., Taylor, K., & Turtle, A. J. 1988, MNRAS, 
234, 155 

Fabbiano, G. 2002, Chemical Enrichment of Intracluster 

and Intergalactic Medium, 253, 93 
Fiore, F., et al. 2001, MNRAS, 327, 771 
Forbes, D. A., & Norris, R. P. 1998, MNRAS, 300, 757 
Galliano, E., Alloin, D., Pantin, E., Lagage, P. O., & Marco, 

O. 2005, A&A, 438, 803 
Galliano, E., Alloin, D., Pantin, E., Granato, G. L., Delva, 

P., Silva, L., Lagage, P. O., & Panuzzo, P. 2008, A&A, 

492, 3 

Heckman, T. M., Dahlem, M., Eales, S. A., Fabbiano, G., 

& Weaver, K. 1996, ApJ, 457, 616 
Hjelm, M., & Lindblad, P. O. 1996, A&A, 305, 727 
Iyomoto, N., Makishima, K., Fukazawa, Y., Tashiro, M., & 

Ishisaki, Y. 1997, PASJ, 49, 425 
Kennicutt, R. C, Jr. 1998, ApJ, 498, 541 
Komossa, S., & Schulz, H. 1998, A&A, 339, 345 
Kristen, H., Jorsater, S., Lindblad, P. O., & Boksenberg, 

A. 1997, A&A, 328, 483 
Lamer, G., Brunner, H., & Staubert, R. 1996, A&A, 311, 

384 

Lindblad, P. O. 1999, A&Ar, 9, 221 

Lonsdale, C. J., & Helou, G. 1985, Cataloged Galaxies and 
Quasars Observed in the IRAS Survey, Jet Propulsion 
Laboratory, Pasadena 

Lutz, D., Maiolino, R., Spoon, H. W. W., & Moorwood, 
A. F. M. 2004, A&A, 418, 465 

Mannucci, F., et al. 2003, A&A, 401, 519 

Massaro, E., Perri, M., Giommi, P., Nesci, R., & Verrecchia, 
F. 2004, A&A, 422, 103 

Massaro, E., et al. 2003, A&A, 399, 33 

Mauch, T., Murphy, T., Buttery, H. J., Curran, J., Hun- 
stead, R. W., Piestrzynski, B., Robertson, J. G., & Sadler, 
E. M. 2003, MNRAS, 342, 1117 

Moorwood, A., et al. 1998, The Messenger, 94, 7 

Ogle, P. M., Marshall, H. L., Lee, J. C, & Canizares, C. R. 
2000, ApJl, 545, L81 

Ogle, P. M., Brookings, T., Canizares, C. R., Lee, J. C, & 
Marshall, H. L. 2003, A&A, 402, 849 

Ott, J., Weifi, A., Henkel, C, & Walter, F. 2005, AIP 
Conf. Proc. 783: The Evolution of Starbursts, 783, 141 

Perlman, E., et al. 2005, ApJ, 625, 727 

Perri, M., et al. 2003, A&A, 407, 453 

Piccinotti, G., Mushotzky, R. F., Boldt, E. A., Holt, S. S., 
Marshall, F. E., Serlemitsos, P. J., & Shafer, R. A. 1982, 
ApJ, 253, 485 

Ranalli, P., Comastri, A., & Setti, G. 2003, A&A, 399, 39 
Read, A. M. 2006, Populations of High Energy Sources in 

Galaxies, 230, 129 
Rector, T. A., Stocke, J. T., Perlman, E. S., Morris, S. L., 

& Gioia, I. M. 2000, AJ, 120, 1626 
Risaliti, G., Elvis, M., Fabbiano, G., Baldi, A., & Zezas, A. 

2005, ApJl, 623, L93 
Risaliti, G., Bianchi, S., Matt, G., Baldi, A., Elvis, M., 

Fabbiano, G., & Zezas, A. 2005, ApJl, 630, L129 



12 



Risaliti, G., Elvis, M., Fabbiano, G., Baldi, A., Zezas, A., 

& Salvati, M. 2007, ApJl, 659, Llll 
Risaliti, G., et al. 2008, MNRAS, L143 
Sakamoto, K., Ho, P. T. P., Mao, R.-Q., Matsushita, S., & 

Peck, A. B. 2007, ApJ, 654, 782 
Sandqvist, A., Joersaeter, S., & Lindblad, P. O. 1995, A&A, 

295, 585 

Sandqvist, A., Elfhag, T., & Jorsater, S. 1988, A&A, 201, 
223 

Schurch, N. J., Warwick, R. S., Griffiths, R. E., & Kahn, 

S. M. 2004, MNRAS, 350, 1 
Schulz, H., Knake, A., & Schmidt-Kaler, T. 1994, A&A, 

288, 425 

Schulz, H., Komossa, S., Schmitz, C, Mucke, A. 1999, 

A&A, 346, 764 
Silich, S., Tenorio-Tagle, G., & Anorve-Zeferino, G. A. 

2005, ApJ, 635, 1116 
Soria, R., Baldi, A., Risaliti, G., Fabbiano, G., King, A., 

La Parola, V., & Zezas, A. 2007, MNRAS, 379, 1313 
Stevens, I. R., Forbes, D. A., & Norris, R. P. 1999, MNRAS, 

306, 479 

Stocke, J. T., et al. 1991, ApJs, 76, 813 
Storchi-Bergmann, T., & Bonatto, C J. 1991, MNRAS, 
250, 138 

Strateva, I. V., Komossa, S. 2009, ApJ, 692, 443, 2009 
Strickland, D. K., Heckman, T. M., Colbert, E. J. M., 

Hoopes, C G., & Weaver, K. A. 2004, ApJ, 606, 829 
Spergel, D. N., et al. 2003, ApJs, 148, 175 
Turner, T. J., Urry, C. M., & Mushotzky, R. F. 1993, ApJ, 

418, 653 

Urry, C. M., et al. 2000, ApJ, 532, 816 

van Buren, D., & Greenhouse, M. A. 1994, ApJ, 431, 640 

Veron, P., Lindblad, P. O., Zuiderwijk, E. J., Veron, M. P., 
& Adam, G. 1980, A&A, 87, 245 

Veilleux, S., Shopbell, P. L., Rupke, D. S., Bland- 
Hawthorn, J., & Cecil, G. 2003, AJ, 126, 2185 

Voges, W., et al. 1999, A&A, 349, 389 

Wang, Q. D., Immler, S., Walterbos, R., Lauroesch, J. T., 
& Breitschwerdt, D. 2001, ApJl, 555, L99 

Wang, Q. D. 2003, High Energy Processes and Phenomena 
in Astrophysics, 214, 32 

Wang, J., Fabbiano, G., Elvis, M., Risaliti, G., Mazzarella, 
J. M., Howell, J. H., & Lord, S. 2009. larXiv:0901.0297l 

Ward, M. J., Done, C, Fabian, A. C, Tennant, A. F., & 
Shafer, R. A. 1988, ApJ, 324, 767 



APPENDIX A: MS 0331.3 3629 

Al Chandra and XMM-Newton Imaging 
Spectroscopy of MS 0331.3 3629 

We present a series of high-quality (up to ~ 2 x 10 4 
counts) X-ray imaging spectra of the high-peaked BL 
Lac (HBL) MS 0331.3-3629. MS 0331.3-3629 (03:33:12.2 
—36: 19:48.0) is a z = 0.308 blazar candidate (|Stocke et al.l 
Il99ll ) which was serendipitously detected in six recent 
long exposures of NGC 1365 - falling within one Chan- 
dra and five XMM-Newton fields. Table IA1I shows a sum- 
mary of the recent Chandra and XMM-Newton observa- 
tions of MS 0331.3— 3629. We used standard Chandra and 
XMM-Newton pipeline reductions to obtain X-ray spectra 



Table Al. MS 0331.3-3629: Chandra/ XMM-Newton Obser- 
vation Summary. (1) Observation ID; (2) Modified Julian Day; 
(3) Total counts; (4) Effective exposure time, after corrections 
for flaring and chip position. 



ObsID 


MJD 


Counts 


T ett 


(1) 


(2) 


(3) 


(4) 






0.3-10 keV 


10 3 s 






ACIS-I2 


ACIS-I2 


3554 


52632 


2758 


14.6 






pn MOS1 MOS2 


pn MOS1 MOS2 


0151370101 


52655 


5888 2767 2700 


12.7 17.4 17.4 


0151370201 


52679 


... 712 675 


... 5.1 5.0 


0151370701 


52864 


... 843 754 


... 8.1 8.2 


0205590301 


53021 


18088 6006 6631 


48.7 57.3 57.4 


0205590401 


53210 


8700 4182 3713 


28.6 46.4 47.1 



(grouped into bins of at least 20 counts per bin) of the 
blazar from an elliptical region with a semi-major axis of 
90" for the pn and 60" for the MOS cameras (for an off 
axis angle ~ 12.5'). About half of the exposure time of Ob- 
sID 0205590401 is lost to flares, while for ObsID 0151370201 
and ObsID 0151370701 the BL Lac falls within the chip-gaps 
of the pn camera, and we use the MOS data only. 

The 0.3-10 keV spectra can be fit by single power-law 
models modified by Galactic absorption. The fit parameters 
are given in the top portion of Table |A"T1 The fits to the two 
highest quality XMM-Newton spectra can be improved by 
adding an intrinsic absorber jV H .i ~ 2 x 10 20 cm~ 2 (F-test 
probability of 1 x 10" 7 and 2 x 10~ 3 for ObsID 0205590301 
and ObsID 0205590401, respectively; last two rows of Ta- 
ble lAlfl . The best spectral fits of all spectra (including intrin- 
sic absorption in the case of the two longest XMM-Newton 
exposures) are shown in Figure IAH Broken power-law fits 
(bknpow in XSPEC) are also consistent with the XMM- 
Newton observations (an F-test shows no improvement to 
the model fit for the Chandra observation) and the fit re- 
sults are presented in Table lA"3l For the five XMM-Newton 
fits, the break energy is around (inbreak) ~ l.lkeVQ the soft 
photon index (ri) ~ 1.9, and hard photon index (T2} ~ 2.1, 
although all values show large variations between the differ- 
ent observations. The current data cannot distinguish be- 
tween the intrinsic absorption and broken power-law fit sce- 
narios; the broken power-law fits have one extra parameter 
but do not provide statistically better fits in comparison to 
the intrinsically-absorbed fits, except for ObsID 0205590401 
(an F-test probability of 2 x 10~ 3 ). 

A2 X-ray Spectral Curvature of MS 0331.3 3629 

iPerlman et al.l (|2005l ) suggest that a large fraction of BL 
Lacs show spectral curvature in the 0.5-10.0 keV band, 
which can be detected using spectral fits in consecutive en- 
ergy bands. This spectral curvature requires a change in 
the standard energy emis sion mechanism (whi ch produces 
a simple power-law), and IPerlman et af, (2005) propose to 
explain this curvature by an episodic (or time-variable) 

7 () indicate weighted averages. 
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Table A2. MS 0331.3-3629: 0.3-10 keV Band Single Power-law Spectroscopic Fits. (1) Obser- 
vation ID; (2) Hard-band rest-frame absorption-corrected ACIS-S or MOS/pn average fluxes in 
units of 10 — 13 ergs -1 cm -2 ; (3) Soft-band rest-frame absorption-corrected ACIS-S or MOS/pn 
average fluxes in units of 10 — 13 ergs -1 cm -2 ; (4) Power-law photon index; (5) Intrinsic abrop- 
tion in units of 10 20 cms — 1 (6) Model fit: "PL" indicates a power law fit [wabs*zpow in 
XSPEC], "APL" indicates a power law fit including intrinsic absorption [wabs*(zwabs*zpow) 
in XSPEC]. A cold absorber equal to the Galactic value, 7V H = 1.4 X 10 20 cm -2 , is included 
in the fits in all cases; (7) Data used in the model fits; (8) \ 2 per degree of freedom (DoF). 



ObsID F 2 -iokcv -Fc>5-2kov T 2V H ,i Model Data x 2 / DoF 

(1) (2) (3) (4) (5) (6) (7) (8) 



3554 


19.6± ;1 




2.06±0.04 


PL 


ACIS-S 


123/111 


0151370101 


20.8±°;7 




1.95±0.02 


PL 


pn+MOS 


490/462 


0151370201 


15.3±j[;| 


".all:* 


2.00±0.05 


PL 


MOS 


60/63 


0151370701 


16.411;? 


17.0+1* 


2.08±0.04 


PL 


MOS 


78/77 


0205590301 


14.9±°o 5 6 


18.11P 


1.97±0.01 


PL 


pn+MOS 


967/953 


0205590401 


i6.2±S:g 


16.6l°-j 


2.08±0.01 


PL 


pn+MOS 


622/664 



0205590301 16.2l°;l 17.6±J;^ 2.05+0.02 2.4+0.5 APL pn+MOS 939/952 
0205590401 17.41°;* 16.0l ° 2.16+0.02 2.0+0.1 APL pn+MOS 612/663 



particle acceleration process. Using the highest signal- 
to-noise best XMM-Newton spectrum of MS 0331.3—3629 
(ObsID 0205590301 ), we looked for signs of curvature simi- 
lar to those found bv lPerlman et all l|2005l) . MS 0331.3-3629 
shows no clear signs of spectral curvature in the 0.5-10 keV 
range. In addition, with the exception of Obs ID 0205590401, 
the broken power-law fit results, which are not preferable 
to a simple absorbed power law, suggest that varying the 
photon index with energy bandpass is not necessary to de - 
scribe the X-ray spectra. Following iPerlman et al.T (2005). 
this could mean that the radio to soft-X-ray emission of 
MS 0331.3— 3629 is consistent with synchrotron radiation 
produced by electrons with a simple power-law distribution. 



A3 Previous X-ray observations of 
MS 0331.3 3629 

MS 0331.3— 3629 was serendipitously detected in eight 
ROSAT observations (see Table lA"4)l . We performed source 
detection and photometry on the 4 PSPC (channels 52-201, 
corresponding to 0.5-2.0 keV) and 3 HRI (channels 0-15, 
0.1-2.4 keV) event files using EXSAS. The resulting count 
rates and 0.5-2.0 keV fluxes (assuming T = 2) are given in 
Table El The RASS detection quoted in Table lA4l was ob- 
taine d from the RO S AT b right source catalog l|Voges et al.l 
1 19991 ). lLamer et all (|l996l ) used the longest PSPC expo- 
sure to measure a photon index in the 0.1-2.4 keV band, 
r soft = 2.1 + 0.15. Table [M] gives the soft-band (0.1-2.5 keV 
in the case of ROSAT PSPC) spectral fits obtained with 
XSPEC in a manner identical to the one used for the Chan- 
dra/ XMM-Newton X-ray spectra. We used standard EXSAS 
processing to obtain the spectra, varying the size of the ex- 
traction region to include 90% of the flux in each case. A 
simple power-law fit (indicated with "PL" in Table |A"5)) with 
absorption equal to the Galactic value and a photon index 
between 1.8 and 2.1 is sufficient to represent the soft-band 
ROSAT PSPC spectra in three of the four cases. In the case 
of RP700921A01, the PL fit is only marginally acceptable 
(x 2 /DoF = 53/40, null hypothesis probability of 8%). The 



inclusion of intrinsic absorption does not improve the model 
fits. 

MS 0331.3-3629 was also observed by BeppoSAX 
(S/N = 3.3) and is pa r t of t he HELLAS sample i|Fiore et al.1 
l200ll V iRector et al] (|2000h show the optical spectrum of 
MS 0331.3— 3629 as part of the Enstetn medium sensitivity 
survey BL Lac catalog. The Enstein flux was fo.3-3.5kov = 
5.22 x 10 -13 ergs -1 cm -2 , consistent with the values listed 
in Table [£J 

A4 X-ray variability of MS 0331.3 3629 

The MS 0331.3— 3629 X-ray flux is variable on timescales 
of weeks to years. The soft (0.5-2.0 keV) band variability 
is about a factor of 2 over 7.5 years (rest-frame); the hard 
(2-10 keV) band variations appear smaller, but the time in- 
terval probed is also shorter. Only the Chandra observa- 
tion shows short term count-rate variability of about 20% 
on timescale of hours, but the evidence for this short term 
variability is not consistent: the Kolmogorov-Smirnov test 
confirms it with a 97% confidence, but a test gives a 
~60% probability of constancy. The power-law photon in- 
dex also varies between the observations on timescales of 
days to years. 

A5 Multiwavelength Observations and the 
Nature of MS 0331.3-3629 

The optical spectrum of MS 0331.3—362 9 shows no emissio n 
lines and has a Ca II H&K break of 30% (|Stocke et al.lll99ll ). 
consistent with a BL Lac classification at 2 = 0.308. An 
HST observation, taken on 1996 October 29, revealed an 
elliptical host which outshines the active nucleus at 7020A 
(the host apparent magnitude is R = 17.83, compa red to 
R = 19.03 for the active nucleus. lUrrv et aL|[2000l ). Note 
that most prior work on the optical-to- X-ray and optical- 
to-radio flux ratios used the total galaxy magnitude in the 
optical, leading to incorrect estimates. The radio-to-X-ray 
index, defined as a rx = - log[F 5Q H Z /-Fikev]/ log[^ 5 GH Z /^ikov], 
is a rx = 0.55. 



14 



MS 0331.3— 3629 has been observed in the J, H, and 
K bands (2 MAS S, likely dominated by the host galaxy, 
Chen et all 120051 ). and the radio: th e VLA (Istocke et al I 



199J >. NVSS (observed on 1993-10, ICondon et all 1 1998 
SUMMS (1998/12-2002/08 composite lMauch et aljhooa i. 
If we construct a spectral energy distribution (SED) based 
on these heterogeneous and non-concurrent observations, 
it peaks at log(i/p) « 17.0, and MS 0331.3- 3629 has the 
appe arance of a high energy peaked blazar l|Donato et al.l 
l200ll ). Assuming that the radio to soft X-ray emission 
is dominated by synchrotron emissio n, a log-parabolic fit 
to th e SED is possible (see for e.g., iMassaro et al.l 120031, 
120041 ). Under the log-parabolic model, the bolometric lu- 
minosity of the s ynchrotron compo n ent can be estimate d 
analytically (e.g., iPerri et alj 120031 ; IMassaro et al. 

I l2004 ). 

L Synch «4 x lO^ergs" 1 . 
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Table A3. MS 0331.3-3629: 0.3-10 kcV Band Broken Power-law Spectral Fits. (1) 
Observation ID; (2) Hard-band rest-frame absorption-corrected ACIS-S or MOS/pn 
average fluxes in units of 10 -13 ergs -1 cm -2 ; (3) Soft-band rest-frame absorption- 
corrected ACIS-S or MOS/pn average fluxes in units of 10 -13 ergs -1 cm -2 ; (4) First 
and (5) Second power-law photon indices; (6) Break energy, in keV; (7) x 2 P er degree 
of freedom (DoF). 



ObsID 


F2-W koV 


-Fo.5-2 koV 




r 2 


-S'brcak 


X 2 /DoF 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


3554 


20.6+"' 2 


15.4+1* 


2.16±0.08 


1.95±0.11 


1.7±0.7 


120/109 


0151370101 


19 3+° :4 


16.0+?; 7 , 


1.91±0.04 


2.01±0.04 


1.3±0.5 


486/460 


0151370201 


16.51 3 ; 2 . 


15-4li; 8 o 


1.87±0.12 


2.2±0.2 


1.4±0.6 


57/61 


0151370701 


16.9+^ 3 


15.5±™ 


2.03±0.06 


2.3±0.3 


2.4±1.5 


76/75 


0205590301 


1 5 2+ '' 9 


17 l +0 - 7 

1, - 1 -0.7 


1.87±0.03 


2.05±0.02 


1.1±0.1 


937/951 


0205590401 


16 6+ ' 9 


15 2 +0 - 7 


1.99±0.03 


2.20±0.04 


1.2±0.2 


605/662 



Table A4. MS 0331.3-3629: ROSAT Photometry. (1) ROSAT observation ID, where 
"RP" refers to the PSPC, "RH" to the HRI, and "RS" to the ROSAT all sky survey; 
(2) Observation date; (3) Effective exposure time; (4) Off-axis angle; (5) Count rate 
in the observed 0.5-2.0 kcV band (PSPC) or 0.1-2.4 keV band (HRI); (6) Rest-frame 
0.5-2.0 keV absorption-corrected flux. 



ObsID 


Date 


T e ff 


OffAngle 


CR 


-Po.5-2 kcV 


(1) 


(2) 


(3) 


(4) 


(5) 






yy-mm-dd 


s 


/ 


counts s 1 


10 13 ergs 1 cm 2 


RP800301N00 


1992-08-22 


840 


43.1 


0.087±0.011 


10±1.3 


RP700921N00 


1992-08-23 


2323 


12.4 


0.090±0.006 


11±0.7 


RP800301A01 


1993-02-02 


4121 


43.1 


0.048±0.004 


5.7±0.5 


RP700921A01 


1993-02-10 


6908 


12.4 


0.062±0.003 


7.5±0.4 


RH701297N00 


1994-08-04 


9182 


12.4 


0.037±0.002 


7.9±0.5 


RH701297A01 


1994-08-23 


323 


12.4 


0.032±0.011 


6.8±2.3 


RH701297A02 


1995-07-04 


9074 


12.4 


0.037±0.002 


7.9±0.5 


RS932309N00 


1990-07-30 


172 




0.090±0.028 


5.2±1.7 



Table A5. MS 0331.3-3629: Soft-band Spectroscopic Fits. (1) Observation ID; (2) Soft- 
band power-law photon index; (3) Rest-frame 2keV monochromatic flux density; (4) Rest- 
frame absorption-corrected fluxes; (5) Bandpass used in the model fit; (6) Model fit: "PL" 
indicates a power law fit [wabs*zpow in XSPEC], "APL" indicates a power law fit including 
intrinsic absorption [wabs*(zwabs*zpow) in XSPEC] with jV H .i = 2.0 ± 0.5 X 10 20 cm -2 . A 
cold absorber equal to the Galactic value, Nu = 1.4 X 10 20 cm -2 is included in the fits in all 
cases; (7) x 2 P er degree of freedom (DoF). 



ObsID 


r so ft 


/2keV 


Fo.S-2 koV 


Data 


Model 


X 2 /DoF 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 






10 -s Jy 


10 13 ergs 1 cm 2 


keV 






RP800301N00 


1.81±0.28 


2.9±0.7 


17+2 


0.1- 


-2.4 


PL 


6/7 


RP700921N00 


1.91±0.09 


2.9±0.3 


18+ 2 


0.1- 


-2.4 


PL 


31/37 


RP800301A01 


1.91±0.12 


1.8±0.3 


111 3 


0.1- 


-2.4 


PL 


19/20 


RP700921A01 


2.10±0.05 


1.9±0.2 


u+j 


0.1- 


-2.4 


PL 


53/40 


3554 


2.12±0.06 


2.7±0.2 


19.9+^ 


0.2- 


-2.5 


PL 


97/87 


0151370101 


1.92±0.02 


2.7±0.1 


15.9+11 


0.2- 


-2.5 


PL 


427/411 


0151370201 


1.87±0.06 


2.5±0.2 


15.ll 2 ; 3 


0.2- 


-2.5 


PL 


59/55 


0151370701 


2.05±0.06 


2.4±0.2 




0.2- 


-2.5 


PL 


72/71 


0205590301 


2.04±0.03 


2.3±0.1 


16.0+i :i i 


0.2- 


-2.5 


APL 


802/784 


0205590301 


1.93±0.01 


2.3±0.1 


u.s+r 6 


0.2- 


-2.5 


PL 


821/785 


0205590401 


2.06±0.02 


2.3±0.1 


16.31S1 


0.2- 


-2.5 


PL 


517/573 
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Figure Al. The 0.3-10 keV spectral fits of MS 0331.3— 3629. In each case, the top panel shows the spectrum and best fit model from 
Table lAll as well as the ObsID, the bottom panel - the fit residuals. 



